Consistency of ground state and spectroscopic measurements on flux qubits 
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We compare the results of ground state and spectroscopic measurements carried out on superconducting flux 
qubits which are effective two-level quantum systems. For a single qubit and for two coupled qubits we show 
excellent agreement between the parameters of the pseudospin Hamiltonian found using both methods. We 
argue, that by making use of the ground state measurements the Hamiltonian of N coupled flux qubits can 
be reconstructed as well at temperatures smaller than the energy level separation. Such a reconstruction of a 
many-qubit Hamiltonian can be useful for future quantum information processing devices. 
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Quantum systems are generally characterized by spectro- 
scopic measurements: the system is excited by electromag- 
netic radiation, with a frequency which matches the level 
spacing, and the response of this excitation is detected. On 
the other hand, quantum theory predicts that the Hamiltonian 
of some quantum-mechanical systems can be completely re- 
constructed from their ground-state properties. For instance, 
quantum mechanical treatment of the ammonia molecule in 
a two-level approximation, shows that its ground state con- 
tains information about time-independent Hamiltonian pa- 
rameters [1]. Superconducting qubits are also described by 
a similar Hamiltonian [2]. They are micrometer- size quantum 
systems [3] which can be easily accessed by a macroscopic 
measuring device. For example, the Hamiltonian parameters 
of a superconducting flux qubit [4] can be determined from 
the measurement of its magnetic susceptibility in the ground 
state (Sfl. In this Letter we will demonstrate that for a single 
and two coupled flux qubits the ground- state and the spectro- 
scopic measurements give the same results. 

The persistent current, or flux, qubit is a small supercon- 
ducting loop with three submicron Josephson junctions [4]. 
Due to the flux quantization only two phases are independent. 
Thus, the circuit is characterized by a two dimensional poten- 
tial U{(j)i, (j)2) which, for suitable qubit parameters, exhibits 
two minima. In the classical case these minima correspond to 
clockwise and anticlockwise supercurrents in the loop. If the 
applied magnetic flux equals half a flux quantum, = ^o/2 
(<l>o = h/2e), both minima have the same potential, leading 
to a degenerate ground state. According to the quantum me- 
chanics the degeneracy is lifted close to this point and the flux 
qubit can be described by the Hamiltonian [2] : 



H(t) = — Ada; — ecFz + Aco^{uot)cFz 



(1) 



where cix , (Jz are the Pauli matrices for the spin basis and 
A is the tunneling amplitude. The qubit bias is given by 
e = Ip{^x — *^o/2), where Ip is the magnitude of the qubit 
persistent current. The last term describes the microwave ir- 
radiation necessary for the spectroscopy which aims to probe 
the stationary energy levels represented by the first two, time 




FIG. 1: Scanning electron micrograph showing two coupled alu- 
minium flux qubits, placed inside a niobium tank coil and do and 
microwave lines. The qubits are fabricated on an oxidized silicon 
substrate by making use of electron-beam lithography and shadow 
evaporation. The loops share a large Josephson junction visible in 
the center. This junction provides a coupling between the qubits with 
a coupling energy J/h = 1.9 GHz. The qubit loops are interrupted 
by three junctions each. Two of them, located in the inner sides, have 
nominal areas of 200 x 700 nm while the outer junctions are 70% 
smaller. The qubit system is placed inside a 30 turn superconduct- 
ing niobium pancake coil which forms a resonant tank circuit, with a 
resonance frequency of 20.8 MHz, with a capacitor mounted on the 
sample holder. The resonator typically has a quality Q = 300. In- 
dividual control of the dc-flux bias of the individual qubits ($a, ^b), 
is provided by the niobium bias lines visible on top of the coil wind- 
ings. The antenna used for supplying the MW-excitation is visible at 
the top of the picture. 



independent, terms of this Hamiltonian. The eigenvalues of 
the Hamiltonian ([T]) depend on the flux bias <l> = <l>a; — $o/2: 



E±($)=±./(/^$)2 + A2. 



(2) 



Spectroscopic measurements detect the excitation of a qubit 
close to the point where the microwave irradiation matches 
the level spacing: hcj = AE{^) = E^{^) - E_{^). By 
measuring AE as a function of $ the qubit parameters A and 
Ip can be obtained as has been shown by van der Wal et al. |0] • 
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Alternatively, the same information can be obtained from 
ground- state measurements. Indeed, let us consider a flux 
qubit weakly coupled to a classical oscillator consisting of an 
inductor Lt and a capacitor Ct forming a tank circuit ITJ. 
Due to the mutual inductance M the tank biases the qubit re- 
sulting in <l> = ^dc + ^rf- Provided that the resonant fre- 
quency of the tank is small, cjt = 1 / ^/LtCt <C A/fi, and 
the temperature is low enough, ksT <C 2 A {ks is Boltz- 
mann's constant), the qubit will reside in its ground state E-. 
The dynamic behavior of the tank-qubit arrangement can be 
described by the Lagrangian 

where q is the charge on the tank capacitor and q is the circu- 
lating current in the tank. Such Lagrangian would lead to the 
nonlinear equation of motion: 

^d^E_{^dc^Mq) 



C = T 



u = ]^LTe 



E_{^dc^Mq) 



^_ ±dC_dC _ 
dt dq dq 



(4) 

however for small amplitude of q the Lagrangian can be lin- 
earized JH] by replacing E_ (g) by its second order Taylor ex- 
pansion around ^dc- Consequently, we obtain the simple La- 
grangian of a particle in a parabolic potential well: 



C= -m q^ 



-'- 7 * 2 



(5) 



The equation of motion which can be obtained from this La- 
grangian is just the simple equation for a particle in a parabolic 
potential, m*q = k*q, where 



M 



(6) 



is the effective mass and k* = 1/ Ct is the curvature of the 
parabolic potential well. Thus, the resonant frequency of the 
tank-qubit arrangement 




— ~ 1 



M2 £E_{<i>ac) 



2Li 



(7) 



contains information on the curvature of the ground state of 
the qubit. Differentiating Eq. Q results in: 



(8) 



(£2($,,)+A2)3/2 ' 

showing that A and Ip can be determined from the depen- 
dence of the resonance frequency of the tank circuit on the 
applied flux. 

In order to compare both methods we fabricated a two qubit 
sample like the one shown in (Fig. 1). As either one of the 
qubits can be biased far away from degeneracy, the single 
qubit properties can be studied as well. This can be under- 
stood if we consider the Hamiltonian of two coupled flux 
qubits: 

(9) 
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FIG. 2: Comparison of the ground state and spectroscopic measure- 
ments for qubit h. (a) Ground state measurements. Presented is the 
dependence of the phase shift between the tank circuit voltage and 
bias current on the flux bias. The solid lines are experimental data 
fitted by the theoretical curves (dashed curves) for qubit parameters 
Ip — 225 nA and A//i = 1.75 GHz. The curves correspond to 
various values of the r/-bias current on the tank circuit resulting in 
rf voltage amplitudes, from top to bottom, Vt — 4.3, 2.9, 0.5, and 
0.3 /iV. (b) Amplitude of the tank voltage as a function of the nor- 
malized magnetic flux in the qubit at the driving frequencies, from 
top to bottom, a;/27r = 18, 5, and 3.5 GHz. The curves have 
been shifted for clarity. The resonant excitation in the flux qubit 
results in the peak-and-dip at the positions defined by the condition 
/\E{<^dc) = (c) Energy gap A^ between the qubit energy 

levels determined from the positions of the mid points of the peak- 
and-dip structures (solid squares). The solid line is the theoretical 
curve calculated from Eq. ^ using the parameters Ip — 225 nA and 
A//i = 1.75 GHz obtained from the ground- state measurements. 
The effective temperature T 70 mK?^ 1.4 GHz-h/ks is smaller 
than the minimal energy level separation 2 A. 



where J is the Josephson coupling energy provided by the 
large connecting Josephson junction. Suppose qubit a is the 
one biased far from its degeneracy point in such a way that Sa 
is large in comparison with the other energy variables. Then, 
qubit a has a well defined ground state with averaged spin 



variables (a 



1 and (a^ 



which can be av- 
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eraged out of the two-qubit Hamiltonian ^ reducing it to: 
H2qbs,red = — A^di^^ — {sb — J)(y^z^ • Apart from the offset in 
the bias term due to the coupHng this is identical to the single 
qubit Hamiltonian ([T]). This offset can be easily compensated 
and measured allowing the determination of the coupling en- 
ergy J [9]. The qubit parameters, A5 and Ip^\ are determined 
from the ground state measurement, as it is described above. 
Analogously, biasing qubit b far from the degeneracy point 
the parameters for qubit a, and /^^^ , can be determined. 
In a similar way the parameters of a A^-qubit Hamiltonian can 
be completely reconstructed from the ground- state measure- 
ments as has already been demonstrated, for instance, for four 
qubits circuits id. 

The qubit parameters can be probed by either the ground 
state (adiabatic) measurements or by making use of spec- 
troscopy. However it naturally raises the question of whether 
the ground-state and the spectroscopy measurements are con- 
sistent? While addressing this problem in this Letter we study 
both approaches in situ for one and two coupled flux qubits. 

Experimentally, the shift of the resonance frequency can be 
obtained by driving the tank circuit with a rf current Irf at 
a frequency close to the resonant frequency ujt and measure 
the phase shift 6 between the rf voltage and driving current. 
For a small qubit inductance L, the phase shift 6 is defined 
byfi]: 



tan O 



M^Q (fE_ 



(10) 



The mutual inductance M, tank inductance Lt and quality 
factor Q can be measured independently giving a value of 
23.4 pH for this pref actor. The results of such measurements 
are shown in Fig. 2(a). Note that the sample was thermally 
anchored to the mixing chamber of a dilution refrigerator at a 
temperature T^ix ^10 mK. The effective temperature of the 
sample T is higher and we estimated from the best theoretical 
fits that T 70 mK[ll]. 

It is important to note that thermal excitations can modify 
the measured signal, which would result in erroneous qubit 
parameters. In practice, thermal excitations are not negligible 
when ksT > 2A. Nevertheless, if ksT < 2A the disper- 
sive measurement provides a correct value of qubit parameters 
[11]. This statement is also confirmed by a good agreement 
between both ground state and spectroscopic measurements 
(see Fig. 2(c)). 

In fact the tank circuit can be used as detector for the spec- 
troscopy measurements as well, since the variation in the pop- 
ulation of the qubits' energy levels results in the change of the 
effective impedance of the tank circuit [12]. The tank circuit 
is insensitive to the microwave signal itself since cot A/fi 
and Q ^ 1. However, if the microwave frequency is close 
to the qubit level separation, the system damps or amplifies 
the voltage on the tank, mimicking the Sisyphus mechanism 
of damping (and heating) of the tank known from quantum 
optics II14I1 . This effect generates the peak-dip structure in 
the Vri^dc) dependence around the resonance (see Fig. 2(b)) 
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FIG. 3: (color online). Landau-Zener interferometry for qubit b. 
Dependence of the tank voltage phase shift G on the dc flux bias 
^dc and the ac flux amplitude $ac (the microwave amplitude). The 
spots along the ^dc axis correspond to the multiphoton resonances 
at the positions defined by the relation AE{^dc) ~ n • hoo; the num- 
bers from 1 to 7 show the position of the n-photon resonances. The 
changes along the $ac axis are due to Stuckelberg oscillations in the 
qubit. The calibration of the driving power of the ac flux can be done 
either with the distance between these oscillations (black arrow) or 
from the slope of the interference fringes (white line). 



12L Il3h . The position of the resonances is the point where a 
peak changes to a dip. From the positions of the mid points 
of the peak-dip structures one can determine the energy gap 
AE between the energy levels. The obtained agreement be- 
tween the adiabatic and spectroscopic measurement for weak 
driving regime is excellent (see Fig. 2(c)). 

With increasing the microwave power, the Landau-Zener 
interference pattern of the qubit is clearly visible. The qubit's 
response in the strong driving regime is demonstrated in Fig. 3 
where the tank voltage phase shift is presented as a function 
of the microwave amplitude and the dc flux bias. The po- 
sition of the multiphoton resonances is approximately given 
by the relation AE{^dc) ^ n ■ hw, where the energy gap 
AE is calculated using the parameters Ip and A obtained 
from the ground state measurement. Moreover, the Landau- 
Zener interferometry allows the calibration of the microwave 
power to the ac flux due to the periodicity of the Stuckelberg 
oscillations on the parameter Alp^ac/^ with the period 27r 



|15|,[16|]. It follows that the distance between the resonances 



(shown by the black arrow in Fig. 3) is approximately equal to 
S^ac = ^Trhuj/Ip. Alternatively, the calibration can be made 
using the slope of the interference fringes (white line in Fig. 3) 
[17, 18]. 

After determining the single qubit parameters far away 
from the degeneracy points, we investigated the two qubit be- 
havior. Firstly, the coupling energy J was determined from 
the offset of the qubit dips from the ^a/b = lines, visible 
in the pure ground state measurements presented in Fig. 4(a). 
Then the qubits were driven by various ac magnetic fluxes 
^ocsincjt. In Fig. 4(b) a frequency in-between both qubit 
gaps was used and therefore only the transitions to the first 
excited state are visible. For higher frequencies, also the sec- 
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FIG. 4: (color online). Spectroscopy of the system of two coupled flux qubits. The dependence of the tank voltage phase shift G on the flux 
biases in qubits a and b is presented for measurements without microwave excitation in (a) and for microwave excitation with a;/27r = 14.125, 
17.625 and 20.75 GHz in (b) till (d) respectively. Inset shows the transition to the third excited level. The blue, magenta and white-dotted 
lines in the pictures with microwave excitation show the expected positions of the resonant excitations of the qubits to the first, second and 
third excited levels respectively, calculated from the energy eigenvalues of Hamiltonian ^ with parameters: Aa/h = 7.9 GHz, Ah/h = 1.75 
GHz, Ip^-* = 120 nA, /p^^ = 225 nA, and J/h = 1.9 GHz. The trough around $b = is due to the ground state curvature of qubit a and 
corresponds to the ground state measurements of Fig. The shallow trough around $a = visible in figure (a), is due to qubit a. 



ond and third excited states become visible as can be seen in 
subfigures (c) and (d). Here also both types of the measure- 
ments (ground-state and spectroscopic) result in the same set 
of parameters for the system. Finally we would like to note 
that the theoretical calculations allow us to plot analogous to 
Figs. 3 and 4 graphs (to be published elsewhere iflll] ). 

In conclusion, the equivalence of the ground-state and spec- 
troscopic approaches for the measurement of the qubit system 
parameters was demonstrated. We have probed the one- and 
two- flux qubit systems by using a dispersive measurement 
technique. It was shown that the ground state measurement 
gives the same qubit parameters as the spectroscopy in the 
weak (Figs. 2 and 4) as well as in the strong driving regime 
(Fig. 3). 
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